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Super-Kamiokande (SK) can search for weakly interacting massive particles (WIMPs) by detecting 
neutrinos produced from WIMP annihilations occurring inside the Sun. In this analysis, we include 
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neutrino events with interaction vertices in the detector in addition to upward-going muons produced 
in the surrounding rock. Compared to the previous result, which used the upward-going muons only, 
the signal acceptances for light (few-GeV/c 2 ~ 200-GeV/c 2 ) WIMPs are significantly increased. We 
fit 3903 days of SK data to search for the contribution of neutrinos from WIMP annihilation in the 
Sun. We found no significant excess over expected atmospheric-neutrino background and the result 
is interpreted in terms of upper limits on WIMP-nucleon elastic scattering cross sections under 
different assumptions about the annihilation channel. We set the current best limits on the spin- 
dependent (SD) WIMP-proton cross section for WIMP masses below 200 GeV/c 2 (at 10 GeV/c 2 , 
1.49xl0~ 39 cm 2 for \X bb and 1.31xlO -40 cm 2 for \X ~t t + t~ annihilation channels), also ruling 
out some fraction of WIMP candidates with spin-independent (SI) coupling in the few-GeV/c 2 mass 
range. 

PACS numbers: 95.35.+d, f4.80.Nb, 14.80.Rt, 96.50.S-, 98.70.Sa 


Weakly interacting massive particles (WIMPs) are fa¬ 
vored as particle candidates for non-baryonic cold dark 
matter (DM), as their interaction strength can explain 
the thermal relic abundance of DM [lH3|. A promising 
way to identify a WIMP DM particle is to search for ex¬ 
cess neutrino flux generated by WIMP self-annihilations 


inside the Sun [4H7] (“WIMP neutrinos”). As the Sun 
travels on the Milky Way arm, WIMPs in the DM halo 
could occasionally become gravitationally bound after 
losing energy by scattering off nuclei in the Sun. The 
WIMPs then pair-annihilate in the deep solar core, pro¬ 
ducing neutrinos from decays of the annihilation prod¬ 
ucts which propagate outward through the Sun and may 
be detected in terrestrial neutrino detectors. An excess of 
a high-energy neutrinos from the Sun, with energies much 
greater than the ~O(10 MeV) solar-fusion neutrinos, ob¬ 
served in neutrino detectors can be interpreted in terms 
of WIMP annihilation rate. We assume that WIMPs 
annihilate via a single channel to a pair of fermions or 
bosons. Also, we assume that the WIMP capture and 
annihilation rates are in equilibrium which enables us 
to remove the WIMP-annihilation cross-section depen¬ 
dence and allows comparison to direct detection results 
on WIMP-nucleon elastic scattering cross section. 

Thanks to the hydrogen-rich composition and large 
gravity of the Sun, tight limits on the spin- 
dependent (SD) scattering cross section of WIMPs on 
protons have been placed by neutrino telescopes such 
as Super-Kamiokande Q, IceCube [f|, Baksan [10j] 
and ANTARES El Recent event excesses or an¬ 
nual modulation signals reported by direct detection ex¬ 
periments such as DAMA/LIBRA [T§|, CoGeNT [lij , 
CRESST 141, CDMS II Si El and conflicting null re¬ 
sults from other direct- and indirect-detection experi¬ 
ments motivate a careful search for light WIMPs be¬ 
low 30 GeV/c 2 . As typical resulting neutrino energies 
are roughly one third to one half of the WIMP mass [1] 
or lower, the Super-Kamiokande (SK) detector’s sensi¬ 
tivity to few-GeV neutrinos makes it suitable for this 

search 53-0. 


Super-Kamiokande is a cylindrical water Cherenkov 
detector located in the Kamioka mine in Japan. The 
inner detector (ID) of 22.5 kton fiducial volume is instru¬ 


mented with 11,129 20-inch Hamamatsu photomultiplier 
tubes (PMTs) and the outer detector (OD) is instru¬ 
mented with 1,885 8-inch PMTs for use as a veto. In¬ 
formation about the experimental setup of the detector 
and its calibration, data reduction and event reconstruc¬ 
tion can be found elsewhere 20|- 22|. 


The SK high-energy neutrino data include neutrino 
events with visible energy>30 MeV and are divided into 
three categories: among contained events in which the 
neutrino interacts inside the ID, fully-contained (FC) 
events have observed Cherenkov light entirely contained 
in the ID, partially-contained (PC) events additionally 
have an exiting particle that deposits energy in the OD. 
Upward-going muons (up-//) are produced by neutrino in¬ 
teractions in the rocks and water surrounding the detec¬ 
tor. FC events are distributed in seven sub-GeV (visible 
energy < 1.33 GeV) and six multi-GeV (visible energy > 
1.33 GeV) sub-categories, and then further divided based 
on particle identification (e-like//x-like), number of re¬ 
constructed Cherenkov rings, number of decay electrons 
and so on. PC events are classified as “OD stopping” 
or “OD through-going” based on their energy deposition 
in the OD. Up-/t’s either stop in the detector (stopping) 
or pass through the detector (through-going). Through- 
going up-// events are further divided into “showering” 
and “non-showering” categories |23| . 

A signal from a 10-GeV/c 2 WIMP will mostly fall in 
the FC sub-GeV and multi-GeV samples, and heavier 
WIMP signals will be mainly distributed in the PC and 
up-/i categories. To increase signal acceptance for light 
WIMPs, FC and PC events are included in the WIMP 
search sample. We take into account all neutrino flavors 
(Ufj,, 17n, i/ e , 17 e , v T and v T ). Compared to the previous 
analysis j§] which used up-// events only, the signal ac¬ 
ceptance has increased 47 times for the 10-GeV/c 2 bb 
channel. 


The analysis in this paper uses data accumulated dur¬ 
ing the SK I-IV run periods. The SK I (1996-2001), II 
(2002-2005, with half PMT coverage), III (2006-2008) pe¬ 
riods correspond to 1489, 799, 518 days of live time for 
the FC/PC samples, and 1646, 828, and 636 days of live 
time for the up-p, sample, respectively. The SK-IV pe¬ 
riod started in 2008; this analysis uses 1096.7 live-days 
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FIG. 1. Angular and reconstructed momentum [GeV/c] distributions of SK I-IV data (black crosses); atmospheric-neutrino 
background MC (normalized to data live time, blue solid); WIMP neutrino signal MC for the 6-GeV/c 2 bb channel (magenta 
dotted) / 200-GeV/c 2 t + t~ channel (cyan dashed) at 90% upper limit, magnified 30 times for visibility. Among seven sub-GeV 
samples, the single-ring e-like 0-decay-electron and /i-like 0,1-decay-electron samples are combined and shown in the left column. 
In the middle column, six multi-GeV samples are combined and shown. PC and up-/r samples are shown in the right column. 


of FC/PC/up-/i data collected until March 2012. 

The Monte Carlo simulation 2^, 241 originally gener¬ 
ated for the primary atmospheric neutrino flux [251 ] is di¬ 
vided into two independent samples and used for two pur¬ 
poses: first, to predict atmospheric-neutrino background, 
and second, to produce the WIMP neutrino signal by 
reweighting to match the signal flux and spectrum. To 
predict the atmospheric-neutrino background, a 250-year 
MC sample per SK run period is normalized to the live 
time of each SK run period and oscillated with parame¬ 
ters: sin 2 9 13 = 0.025, sin 2 d 12 = 0.304, sin 2 9 23 = 0.425, 
A m 21 = 7.66 x 10 -5 eV 2 and Am 32 = 2.66 x 10 -3 eV 2 . 
In the other 250 years of MC sample per run period, 
events coming approximately from the solar direction are 
selected and then weighted by the ratio of WIMP neu¬ 
trino flux to original atmospheric neutrino flux. To sim¬ 
ulate the neutrino flux from WIMP annihilation in the 
Sun at SK, we used the WIMP MC simulator WimpSim 
3.01 J2(|. Focusing on the light WIMP search, we con¬ 
sidered WIMP masses from 4 GeV/c 2 to 200 GeV/c 2 for 
XX -A t + t~ /bb/W + W~ channels. The spectrum of the 
bb channel is “softest” and the t + t~ channel produces a 
“harder” spectrum. 


In order to discriminate the large background of atmo¬ 
spheric neutrinos at low energies, which have an oc E~ 2 * * ' 7 
power-law spectrum, and to take into account the energy- 
dependent angular correlation for low-energy events, we 
perform a least-squares fit which makes full use of the en¬ 
ergy, angle and flavor information. Under the hypothesis 
that the collected SK data consist of both atmospheric 
neutrinos and WIMP neutrinos, we compare data to MC. 
The data and MC in the 18 sub-categories are distributed 
into 480 bins based on reconstructed momentum and 
cos d sun , where 9 sun is the angle between the Sun and the 
reconstructed event direction. Figure[[|shows an example 
of data distributions compared to the expected WIMP- 
induced signal. We fit the data using the least-squares 
technique based on Poisson statistics to find the amount 
of signal contribution added to background which best 
matches the data. For each mass and annihilation chan¬ 
nel, the energy spectrum and flavor composition of signal 
neutrino flux are fixed and the global normalization of the 
signal flux is allowed to vary freely. The “pull” approach 
allows us to incorporate systematic uncertainties [27] in 
the fitting as: 


X = mm 
PM 


2 g { N exp _ N data + jydata^jydata/jyexp^ + £ 

1=1 1=1 

where = N? G (1 + £ ffa) + /3iVf (1 + £ #l), 

l=i l=i 


r 


and i is the index of the bins; Nf ata is the number of 


events observed in each bin; N^ G is the background ex- 
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pectation in the bin; Nf is the number of events of signal 
MC in the bin, where f3, the global normalization param¬ 
eter for signal, stands for the allowed fraction of signal; 
aj is the la value of the j -th systematic uncertainty; tj 
is its “pull” and fj/aj is the predicted fractional 

change of the number of background (signal) events in 
the *-th bin due to a la change of the j-th systematic 
uncertainty. In the pull approach, y 2 is estimated by 

solving Sx 2 /Sej = 0 for all e/s [53 ]. 

For background (signal), 66 (48) sources of system¬ 
atic uncertainties are considered. A total of 16 system¬ 
atic uncertainties related to neutrino interaction in SK 
are considered in common for both background and sig¬ 
nal. Another 14 uncertainty sources related to event re¬ 
duction, and 11 sources related to event reconstruction 
and selection are calculated for SK I, II, III and IV in¬ 
dependently and considered for both signal and back¬ 
ground. For background MC, a total of 25 uncertainties 
related to prediction of atmospheric neutrino flux, includ¬ 
ing la uncertainties on oscillation parameters, are consid¬ 
ered. All the uncertainties in neutrino interaction, event 
reduction, reconstruction, selection and background at¬ 
mospheric neutrino flux are similar to those described 
in [28( except the ones related to zenith angle distribu¬ 
tions. These uncertainties are listed elsewhere (for ex¬ 
ample, Tab. VI,VII,VIII in J28|). For signal MC, 7 un¬ 
certainties related to the oscillation of WIMP neutrinos 
during propagation through the Sun, vacuum, and the 
Earth, are considered. The full systematic uncertainty 
has up to 10% effect on the WIMP sensitivity, where the 
largest contributions come from neutrino interaction and 
atmospheric neutrino flux uncertainties. 



FIG. 2. The 90% upper limit on total integrated muon- 
neutrino flux from WIMP annihilations in the Sun at SK for 
the t + t~ channel in red solid; bb in green dashed; W + W~ 
in blue dot-dashed. The shadowed regions show la bands of 
the sensitivity study results (color scheme is the same as for 
data). 


The best-fit value of /?, f3 m in, is defined as the value 
at which y 2 is minimized with respect to (3. For all 
tested WIMP hypotheses, we found the resulting /3 mirl 
values are all consistent with the hypothesis of no WIMP 
neutrino contribution, and we set 90% upper limits as¬ 
suming that the obtained y 2 values approximately follow 
a normal distribution. To set a physically-meaningful 
confidence limit, we follow a Bayesian approach [29] and 
renormalise the f3 distribution in the physically-allowed 
region (/3 > 0). The fitted y 2 values at minimum (y^ in ; 
DOF is 479), the Ay 2 value at (3 = 0 from the minimum, 
and the Ay 2 values for 90% upper limit (Ayg 0 ) calcu¬ 
lated with the Bayesian approach are shown in Tab. [T] for 
the tested WIMP mass and annihilation channels. 

Figure [2] shows the derived 90% upper limit on the 
muon-neutrino flux from WIMP annihilations in the Sun 
at SK. The limit shown here for is an example; limits 
on other flavors or antineutrinos are determined by ap¬ 
propriate flavor scaling for the given WIMP annihilation 
channel assumption. The sensitivity study was done by 
substituting SK data with 500 sets of toy MC generated 
for the hypothesis of no WIMP-annihilation contribution 
and the result is shown together in Fig. [2] as a band rep¬ 
resenting a ler range of results. 



FIG. 3. 90% C. L. upper limits on SD WIMP-proton cross 
section calculated at DarkSUSY j3o| default are shown in red 
solid with uncertainty bands to take account uncertainties in 
the capture rate for the bb, W + W~ and t + t~ channels from 
top to the bottom. Also shown are limits from other exper¬ 
iments: IceCube Q in brown dashed: bb (top) / W + W~ 
or t + t~ (bottom); BAKSAN 0 in pink dot-dashed: bb 
(top) / W + W~ (middle) / t + t~ (bottom); PICASSO [U 
(blue long-dashed); SIMPLE [33] (green long dot-dashed). 
The black shaded region is the 3a C.L. signal claimed by 
DAMA/LIBRA [Ta, 33]. 

Using DarkSUSY 5.0.6 [30], we convert the upper limit 
on the neutrino flux to upper limits on WIMP-nucleon 
cross sections. We assume that WIMPs have only a sin¬ 
gle type of interaction with a nucleus, either an axial 
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6 
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4.9 
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1.25 

4.9 

70.8 
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44.5 
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10 

bb 
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20 
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50 
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3.73 

100 

bb 
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32.7 
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0.14 
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1.24 

2.76 

1.26 

200 

bb 
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44.9 
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f 

is 
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0.07 
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0.17 

3.00 

4.23 

2.03 


T^T 

508.9 

0.08 

2.4 

0.19 

1.33 

1.88 
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TABLE I. Xmin i Ax 2 at /3 = 0, Axlo (Ax 2 for 90% Bayesian upper limit), 90% upper limit on the muon-neutrino flux 
from WIMP annihilations in the Sun at SK and SD/SI/SI (IVDM) scattering cross section limits for each WIMP mass and 
annihilation channel. 



and protons (/ p ). We consider two representative ex¬ 
amples: the commonly-assumed isospin-invariant case 
{fn/fp = 1) and the case of isospin-violating dark mat¬ 
ter (IVDM) [39j] with destructive interference f n /f p = 


—0.7 [40M431 ] . The latter is calculated based on [4 4 
Also assumed are a standard DM halo with local density 
0.3 GeV/cm 3 [H,|46|, a Maxwellian velocity distribution 
with an RMS velocity of 270 km/s and a solar rotation 
speed of 220 km/s. The results are listed in Tab. Q] They 
are also plotted together with other experimental results 
in Fig. [3] for SD coupling and Fig. 0] for SI coupling for 
the isospin-invariant case. For IVDM, we note that for 
the result given in the Tab.Q] the entire CDMS II Si 151 
90% C.L. signal region will be in tension with our t + t~ 
channel result. 


FIG. 4. 90% C. L. upper limits on the SI WIMP-nucleon cross 
section (plotting scheme is the same as FigO • Also shown are 
event excesses or annual modulation signals reported by other 
experiments: DAMA/LIBRA (black shaded regions, 3<r C.L.); 
CoGeNT [l3j] (ma gen ta diagonally cross-hatched region, 90% 
C.L.); CRESSTII 111 (violet horizontally-shaded regions, 2cr 
C.L.); CDMS II STM (blue vertically-shaded region, 90% 
C.L.); and limits: IceCube jtij] in brown dashed: bb (top) / 
W + W~ or t + t~ (bottom); SuperCDMS 01 (cyan dotted); 
CDMSlite [U] (blue long dot-dashed); XENON10 S2-only 0 
(dark green dash triple dot); XENON100 0] (green dash 
double dot); LUX [0] (orange long-dashed). 


vector interaction in which WIMPs couple to the nuclear 
spin (SD) or a scalar interaction in which WIMPs cou¬ 
ple to the nucleus mass (spin-independent, SI). The SI 
coupling can have different couplings to neutrons ( f n ) 


There are several sources of uncertainties related to the 
WIMP capture process. Uncertainty in the composition 
of the Sun is considered by comparing the DarkSUSY de¬ 
fault choice BS2005-OP model with the BS2005-AGS,OP 
model [47| with lower heavy element abundances. We 
compared the Helm-Gould nuclear form factor [?| used in 
DarkSUSY to the choices in 48| and J49|. These effects 
are predicted to be small in the SD-coupling case and 
up to 25/45% in the Si-coupling where heavier elements 
than hydrogen contribute to the capture. The effect of 
uncertainties in the velocity distribution of WIMPs is 
determined in [50| to be up to 40(25)% for SD(SI) cou¬ 
plings. The effects of the planets on the capture rate 
is determined to be negligible [U, [52:] . Solar evapora¬ 
tion is expected to have no impact above a WIMP mass 
of 4 GeV/c 2 H, 53)-[55|. These uncertainties are added 


in quadrature and indicated by the shadowed regions in 
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Fig. [3] and Fig. [4] Uncertainty in the local WIMP den¬ 
sity will make a similar vertical shift of the limits for all 
direct and indirect detection experiments, and so is not 
indicated. 

In conclusion, the result of the first WIMP search using 
contained events in SK is presented. No significant signal 
excess was found for 4-200-GeV/c 2 WIMP hypotheses. 
The derived upper limit on the SD WIMP-proton cross 
section places the most stringent constraint to date for 
WIMP masses below 200 GeV/c 2 even for the softest 
( bb ) channel, assuming the equilibrium condition between 
capture and annihilation rates. For the Si-coupling case 
with 100% annihilation to r + r _ , we exclude new regions 
for WIMP masses below 6 GeV/c 2 . 
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